y  fl0-fli49  845  SIDE  MOMENT  EXERTED  BV  ft  TUO-COMPONENT  LIQUID  PAVLOAD  V 
ON  A  SPINNING  PROJECTILECU)  flRMV  BALLISTIC  RESEARCH  LAB 
ABERDEEN  PROVING  GROUND  MD  C  H  MURPHV  DEC  84 
UNCLASSIFIED  BRL-TR-2624  SBI-AD-F200  556  F/G  19/1  ML 


MB  file  copy 


AD-A149  845 


B 

R 

L 


TECHNICAL  REPORT  BRL-TR-2624 


SIDE  MOMENT  EXERTED  BY  A  TWO-COMPONENT 
LIQUID  PAYLOAD  ON  A  SPINNING  PROJECTILE 


Charles  H.  Murphy 
December  1984 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


US  ARMY  BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND.  MARYLAND 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies 
from  the  National 
U.  S.  Department 
22161 . 


of  this  report  may  be  obtained 
Technical  Information  Service, 
of  Commerce,  Springfield,  Virginia 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other 
authorized  documents. 


The  use  of  trade  names  or  manufacturers'  names  in  this  report 
does  not  constitute  indorsement  of  any  commercial  product. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  This  PAGE  'When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  REPORT  NUMBER  ,2  GOVT  ACCESSION  NO. 

HSaaDSIBZSiBiBEEEQXHHH 

TECHNICAL  REPORT  BRL-TR-2624  /}/*//■;  c 

4  TITLE  (and  Subtle!#) 

TYPE  of  REpDRT  4  PERIOD  COVERED 

SIDE  MOMENT  EXERTED  BY  A  TWO-COMPONENT 

Final 

LIQUID  PAYLOAD  ON  A  SPINNING  PROJECTILE 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  author(»; 

8  CONTRACT  OR  GRANT  NUMBER.*' 

Charles  H.  Murphy 

9-  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

10  PROGRAM  ELEMENT.  PROJECT,  TASK 

US  Army  Ballistic  Research  Laboratory 

AREA  1  WORK  UNIT  NUMBERS 

ATTN:  AMXBR-LFD 

Aberdeen  Proving  Ground,  Maryland  21005-6066 

RDT&E  1L161 102AH43 

M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

12.  REPORT  DATE 

US  Army  Ballistic  Research  Laboratory 

December  1984 

ATTN:  AMXBR-OD-ST 

13  NUMBER  OF  PAGE > 

79 

14  MONITORING  AGENCY  NAME  4  ADDRES§(7f  ditferoru  from  Controlling  Office) 

15  SECURITY  CLASS  (of  thle  report) 

UNCLASSIFIED 

1S«  DECLASSIFICATION  DOWNGRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  t hie  Report ) 

Approved  for  public  release;  distribution  unlimited 

• 

17.  DISTRIBUTION  STATEMENT  (of  tho  abolrocf  onlorod  In  Block  30,  II  dllloront  from  Report) 

10.  SUPPLEMENTARY  NOTES 

» 

19  KEY  WORDS  (Continue  on  reveree  eide  if  neceeeery  end  Identify  by  block  number) 

Eigenfrequency 

Gyroscope 

Liquid  Moment 

Spinning  Projectile 

Two-Liquid  Payload 

20  ABSTRACT  rCootfoue  on  rowereo  •  if  nmcmeaety  mtd  Identify  by  block  number) 

(bja) 

A  linear  boundary  layer  theory  is  derived  for 

two  liquids  in  a  cylindrical 

cavity  in  a  coning  and  spinning  projectile.  Predicted  eigenfrequencies  agree 

well  with  available  experimental  results.  Side  moment  coefficients  are  com- 

puted  for  various  cases.  It  is  shown  that  significant  increases  in  the  single- 

liquid  side  moment  can  be  caused  by  the  addition  of 

a  small  amount  of  heavier  1 

liquid  and  for  certain  locations  of  the  liquid  interface,  pairs  of  eigenfre-  t 

[quencies  can  coalesce. 

1 

00  /«*,  M73 


EDITION  or  I  NOV  88  IS  OUSOLE  TC 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  or  This  P*  .r  '■»>•*>  I'olo  Fnlrr*./) 


TABLE  OF  CONTENTS 


LIST  OF  ILLUSTRATIONS .  5 


LIST  OF  TABLES. 


I.  INTRODUCTION. 


II.  LIQUID  BOUNDARY  CONDITIONS .  10 

III.  THE  INVISCID  SOLUTION .  14 

IV.  THE  VISCOUS  SOLUTION . . .  16 

V.  INVISCID  BOUNDARY  CONDITIONS .  21 

VI.  LIQUID  MOMENT .  23 

VII.  DISCUSSION .  25 

VIII.  SUMMARY .  27 

ACKNOWLEDGMENT .  28 

REFERENCES .  39 

APPENDIX  A.  EFFECT  OF  CENTRAL  ROD .  41 

APPENDIX  B.  EFFECT  OF  DIFFERENT  KINEMATIC  VISCOSITIES .  49 

APPENDIX  C.  AXISYMMETRIC  EIGENVALUES .  61 

LIST  OF  SYMBOLS .  67 

DISTRIBUTION  LIST .  73 


DTIC 

I^ELECTE 
k  FEB  4  1985 


Accession  For 


NT  IS  OP.’J-.i 
DTIC  T\ 
Unarm  -  c  - 

By— 

JHaV' ' 

Avail'1 


1 


LIST  OF  ILLUSTRATIONS 


versus  b^/a  for  Re^  =  4  x  104,  c/a  =  3.1, 
f  =  1  and  Various  Density  Ratios . 

versus  b^/a  for  Re^  =  4  x  104,  c/a  =  3.1,  f  =  .98 
with  a  Free  Surface  and  Various  Density  Ratios . 

versus  bj/a  for  Re^  =  4  x  104,  c/a  =  3.1,  f  =  .98 
with  a  Central  Rod  and  Various  Density  Ratios . 

versus  bj/a  for  Re^  =  2  x  106,  c/a  =  3.127,  =  .82 

and  Various  Free  Surface  Fill  Ratios.  Experimental  Data  are 
from  Reference  9 . 

t-ji  Maximum  Side  Moment  Coefficient  versus  bj/a  for 
Rej  =  4  x  104,  c/a  =  3.1,  f  =  1,  and  Various 
Density  Ratios . 

t-ji  and  T52  versus  b^/a  for  Re^  =  4  x  104,  c/a  =  3.1, 

f  =  1*  P21  =  . . . . . . . . 

T72  versus  bj/a  for  Re^  =  106,  c/a  =  4.29,  f  =  1 

and  Various  Density  Ratios . 

T72  and  j  versus  bj/a  for  Re^  =  106,  c/a  =  4.29, 

f  =  1  and  P2^  =  0.8 . 

Xj2  Maximum  Side  Moment  Coefficient  versus  b^/a  for 

Re^  =  106,  c/a  =  4.29,  f  =  1  and  Various  Density  Ratios . 

Damping  Rate,  -et,  versus  r  for  Re,  =  5.34  x  105, 


LIST  OF  ILLUSTRATIONS  (cont'd) 


31  Maximum  Side  Moment  Coefficient  versus  b./a  for 

4  * 

Re^  =  4  x  10  ,  c/a=3.1,  f = . 98  with  a  Central  Rod  and 

Various  Density  Ratios . 

,,  Maximum  Side  Moment  Coefficient  versus  b./a  for 
Re^=4  x  10  ,  c/a=3.1,  f=.98,  ^21  =  ‘^  anc*  Various 
Kinematic  Viscosity  Ratios . 

2jq  versus  b^/a  for  Re^=4.3  x  104,  c/a=0.995,  f=l, 
m=0  and  Various  Density  Ratios . 66 


LIST  OF  TABLES 


Page 

Inviscid  Perturbation  Functions.. .  ^ 

Viscous  Coefficient  Functions  for  Internal  Free  Surface .  19 

Viscous  Coefficient  Functions  for  Central  Rod .  44 

_  54 

Pressure  and  Radial  Velocity  Functions . 

Coefficients  of  Eq.  (Bl) .  88 

Equations  for  (EjK,  F^,  E£K,  .  ^ 

Viscous  Radial  Velocity  Functions  for  Central  Rod .  58 

Inviscid  Perturbation  Functions  for  m  *  1 .  64 


7 


I.  INTRODUCTION 


The  prediction  of  the  complete  moment  exerted  by  a  spinning  liquid  pay- 
load  on  a  spinning  and  coning  projectile  has  been  a  problem  of  considerable 
interest  to  the  Army  for  some  time.  For  a  fully  spinning  liquid,  the  linear 
side  moment  was  first  computed  by  Stewartson1  for  an  inviscid  payload  by  use 
of  eigenfrequencies  determined  by  the  fineness  ratio  of  the  cylindrical  con¬ 
tainer.  Wedemeyer2  introduced  boundary  layers  on  the  walls  of  the  container 
and  was  able  to  determine  viscous  corrections  for  Stewartson's  eigenfrequen¬ 
cies,  which  could  then  be  used  in  Stewartson's  side  moment  calculation. 
Murphy3  then  completed  the  linear  boundary  layer  theory  by  including  all 
pressure  and  wall  shear  contri butions  to  the  liquid-induced  side  moment.  The 
Stewartson-Wedemeyer  eigenvalue  calculations  have  been  improved  for  low 
Reynolds  numbers  by  Kitchens  et  at .4  through  the  replacement  of  the  cylindrical 
wall  boundary  approximation  by  a  linearized  Navier-Stokes  approach.  Next, 
Gerber  et  aZ.5”6  extended  this  linearized  NS  technique  to  compute  better  side 
moment  coefficients  for  Reynolds  numbers  less  than  10,000.  Finally,  the  roll 
moment  for  a  fully  spun-up  liquid  was  computed  by  Murphy.7"8 


1 .  K.  Stewartson ,  "On  the  Stability  of  a  Spinning  Top  Containing  Liquid," 
Journal  of  Fluid  Mechanics ,  Vol.  5,  Part  4,  September  1959,  pp.  577-592. 

2.  E.  H.  Wedemeyer,  "Viscous  Correction  to  Stewartson' s  Stability  Criter¬ 
ion,"  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground,  Maryland, 
BRL  Report  No.  1325,  June  1966.  (AD  489687) 

3.  C.  H.  Murphy,  "Angular  Motion  of  a  Spinning  Projectile  With  a  Viscous 

Liquid  Payload,"  Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground, 
Maryland,  BRL  Memorandum  Report  ARBRL-MR-03194 ,  August  1982.  (AD 
A118676 ) .  (See  also  Journal  of  Guidance,  Control,  and  Dynamics,  Vol.  6, 
July- August  1983 ,  pp.  280-286 . ) 

4.  C.  W.  Kitchens,  Jr.,  N.  Gerber,  and  R.  Sedney,  "Oscillations  of  a  Liquid 

in  a  Rotating  Cylinder:  Solid  Body  Rotation ,"  Ballistic  Research  Labora¬ 
tory,  Aberdeen  Proving  Ground,  Maryland,  BRL  Technical  Report  BRL-TR- 

02081,  June  1978.  (AD  A057759) 

5.  N .  Gerber,  R.  Sedney ,  and  J.  M.  Bartos,  "Pressure  Moment  on  a  Liquid- 

Filled  Projectile:  Solid  Body  Rotation,"  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  Maryland,  BRL  Technical  Report  ARBRL-TR-02422 , 
October  1982.  (AD  A120567) 

6.  N.  Gerber  and  R.  Sedney,  "Moment  on  a  Liquid-Filled  Spinning  and  Nutating 
Projectile:  Solid  Body  Rotation ,"  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland,  BRL  Technical  Report  ARBRL-TR-0247 0 ,  February 
1983.  (AD  A125332) 

7.  C.  H.  Murphy,  "Liquid  FxjIj-^I  Do 1 1  Moment  Induced  by  a  Spinning  and 

Coning  Projectile,'  Ballistic  Research  Laboratory,  Aberdeen  Proving 
Ground,  Maryland,  BRL  Technical  Report  ARBRL-TR-0252 1 ,  September  1983. 
(AD  A133681)  (See  also  AIAA  Paper  83-2142,  August  1983.) 

8.  C.  H.  Murphy,  "A  Relationship  Between  Liquid  Roll  Moment  and  Liquid  Side 

Moment,"  Ballistic  Research  Labor itory ,  Aberdeen  Proving  Ground, 

Maryland,  BRL  Memorandum  Report  ARBRL-MR-03347 ,  April  1984.  (AD  A140658) 


An  important  limitation  of  this  work  has  been  the  restriction  to  a  single 
liquid  in  the  payload.  Scott9  has  derived  relations  for  eigenfrequencies  for 
two  inviscid  liquids  and  obtained  fair  experimental  agreement.  For  Scott's 
inviscid  liquids  the  tangential  perturbation  velocities  have  discontinuous 
jumps  at  the  two-liquid  interface. 

In  this  report  we  will  consider  two  viscous  liquids,  although  we  will 
restrict  our  consideration  of  viscosity  to  boundary  layers  near  the  cylin¬ 
drical  walls  and  near  the  two-liquid  interface.  The  remainder  of  the  liquid 
will  be  considered  to  be  inviscid.  Under  these  assumptions,  side  moment 
coefficients, as  well  as  eigenf requencies, wi 1 1  be  computed. 


II.  LIQUID  BOUNDARY  CONDITIONS 

Two  coordinate  systems  will  be  used  in  this  report:  the  nonrolling  aero- 

ballistic  XYZ  system  whose  X-axis  is  fixed  along  the  missile's  axis  of  sym¬ 
metry  and  the  inertial  XYZ  system  whose  X-axis  is  tangent  to  the  trajectory  at 
time  zero.  Both  coordinate  systems  have  origins  at  the  center  of  the  cylin¬ 
drical  payload  cavity,  which  is  assumed  to  be  at  the  center  of  mass  of  the 
projectile.  Location  in  the  cavity  can  be  specified  in  the  aerobal 1 i Stic 

system  by  the  cylindrical  coordinates  x,  r,  e  and  in  the  inertial  system  by  x, 

r,  6.  The  boundary  of  the  cavity  is  given  by  x  =  ±c  and  r  =  a  where  2c  is  the 
height  of  the  cavity  and  2 a  is  its  diameter.  The  projectile  is  assumed  to  be 

performing  a  coning  motion  of  amplitude  K^(t)  and  phase  angle  <^(t).  If  a  and 

3  are  the  angles  of  attack  and  side-slip  of  the  missile's  axis  with  respect  to 
the  trajectory  (the  X-axis), 

-  -  i  4>i 

B  +  i  a  =  Kje  (2.r 


where  <f>  =  <Jt 


s  =  ( e  +  i )  t 

^  i 

K  =  K^O)  e 


9.  W.  E.  Scott,  ''The  Inertial  Wave  Frequency  Spectrum  in  a  Cylindrically 
Confined,  Inviscid,  Incompressible  Two  Component  Liquid,  Ballistic 
Research  laboratory ,  Aberieen  Proviny  Ground,  Maryland,  BRL  Report  No. 
1609,  September  19?2.  (AD  7.5 2439).  (See  also  Physics  of  Fluids,  Vol. 
16,  No.  1,  po.  .9-12,  January  19?3.)  —  -  ~ 


and  $  is  the  axial  component  of  the  angular  velocity  -'elative  to  inertia  axes 
and  is  assumed  to  be  positive  and  constant.  Linear  relations  Detween  cylin¬ 
drical  coordinates  in  the  two  coordinate  systems  were  derived  in  Reference  3: 


x  =  x  -  r  cos 

(  $j-6)  =  x-r  R  'K  eS  ;  f 

(2.2) 

r  =  r  +  x  K^cos 

( $,  - 9)  =  r  +  x  R  -,K  e3  *■  / 

(2.3) 

9=6+  (x/r)  K, 

sin  (;i  -  j)  =  9  -  (x/r)  R  { i K  eSi!}  , 

(2.4) 

where  R{  }  =  [{  }  +  f}]/ 2  is  the  real  part  of  a  complex  quantity. 

The  two  liquids  have  densities  ^  and  ^(pj  >  p£)  and  kinematic  vis¬ 
cosities  and  v^.  When  the  liquids  are  fully  spun  up  and  Kj  =  0,  liquid  1 
occupies  the  cylindrical  annular  region  b^  <  r  <  a  and  liquid  2  occupies  a 
cylindrical  annular  region  b^  <  r  <  bj .  The  fill  ratio,  f,  is  the  ratio  of 
the  volume  of  the  annular  region  containing  both  liquids  to  that  of  the  com- 

p 

plete  cylinder,  i.e.,  1  -  (b^/a)  .  The  surface  r  =  b2  is  either  a  free 

surface  or  the  surface  of  a  rigid  central  rod.  The  free  surface  boundary  will 
be  considered  in  this  report  while  the  central  rod  will  be  discussed  in 
Appendix  A.  For  a  fully  filled  cylinder,  b£  =  0  and  f  =  1. 

When  the  cylinder  is  forced  to  perform  a  coning  motion,  the  interface 
between  the  liquids  is  located  at 

r !  =  bL  (1  +  Ml)  (2.5) 

where 

|  !  , 

and  the  free  inner  surface  is  locatt  1  at 


The  velocity  components  in  the  two  liquids  have  a  very  simple  form. 


v  .  =  R  i.u  .  eS '  1  •  a 
M  js 


V  .  =  R  (v.  eS',-1J}  a. 

hj  JS 


V  .  =  •  r  +  R  {w.  eScrl-J}  a' 
‘J  J  s 


where  j  =  1,2. 

The  fluid  pressure,  however,  has  a  more  complicated  form: 


p  =  p 


0 


0  <  r  <  r- 


Pt)  +  2  * 


2  k  2 

r  -  b2 


R  ip2ses^ie}  .:,2a2;2 


r2  -  r  -  rl 


P0  +  °2  ^ 


.  2  .  1 
B1  ‘  b2 

+  > 

f2 

r  -1 

?  2 

CM 

_ 1 

2  *■ 

- j 

CM 

_ _ J 

♦ R  (pis n  ^ 


r1  <  r  <  a 


(2.71 

(2.8) 

(2.9) 


2.19) 


If  we  assume  the  small  perturbation  of  the  interface  surface  has  the  same  form 
as  the  perturbation  velocities. 


nj  =  R  (n1$(x)  eS  r-i  1  • 


(2.11) 


■' v  ise  of  Eqs.  (2.3,  2.11),  the  equation  of  the  interface  surface  is 
F  (  x , r ,  ! , t )  -  r  -  b,  -  R  fb,  n.  e5*"1’’} 


At  the  interface,  the  two  liquids  have  the  same  velocity,  i.e.,  the 
velocity  of  the  interface  itself. 


dP  }P  ,,  ,-f  ^3)  if  iF  _ 

=  Vxj  17  +  Vj  17  +  ~~  T)  +  W  =  0 


(2.13) 


When  only  linear  terms  are  retained,  this  reduces  to 


a  vls  (blf  x)  a  v?c  (b,  ,x 


2s  ^U1 


bl  nis(x^ 

s  -  i  s  -  i 

Similarly,  the  perturbation  of  the  inner  surface  has  the  form 

a  v2$  (b2,x) 


b2  n2s(x) 


s  -  i 


where 


=  R  (n2s  e 


s  4>-i  o 


(2.14) 


(2.15) 


In  addition  to  three  continuous  velocity  components  at  the  interface,  we 
will  require  the  pressure  and  the  two  viscous  shears  to  be  continuous  at  the 
interface.  Equations  (2.10,  2.14)  yield  a  simple  pressure  relation  at  the 
interface : 


pls 


(bpx)  + 


bl  vis(bi’x) 

a(s  -  i) 


°2 1 


P2s(brx)  + 


bl  v2slbl,x^ 
a  (s  -  i) 


(2.16) 


where  ^21  ”  p 2/ p^. 


If  only  velocity  gradients  normal  to  the  interface  are  considered,  the 
continuous  viscous  shear  assumptions  imply  that 


UlS  (b  1  ,X  )  u?clbl»x) 


Ur 


2sv”r 


(2.17) 


a  r 


3  r 


J  W|,  (  b  p  X  ; 


•  w v . fb, , > ) 


i  r 


>s  _i_ 
r 


\  2  .  i  H  j 


At  the  free  surface,  tin-  press  -re  i  >  constant,  and  l.’u-  •  i  scons  snears  are 


zero. 


(2.19) 


P2s^2,x)  ^ 


b2v2s(b2*x) 
a  (s  -  i) 

3  u2s(b2,x) 


(2.20) 


3  w2$(b2,x) 


(2.21) 


Finally,  at  the  cylinder  surface  the  liquid  must  have  the  rigid  body  motion  of 
the  cylinder. 


u1c  =  (s  -  i)  K 


(2.22) 


'ls  =  -(s  -  i)  (x/a)  K 


(2.23) 


w,  =  i  (s  -  i)  (x/a)  K. 


(2.24) 


III.  THE  INVISCID  SOLUTION 

The  perturbation  functions  are  written  as  sums  of  viscous  and  inviscid 
parts. 


ujs  uj  si  +  ujsv 


wjs  =  wjsi  +  wjsv 


(3.1) 


V  ■  =  V  •  ■  +  V  • 

js  JS1  jsv 


pjs  pjsi  +  pjsv  . 


The  differential  equations  for  the  inviscid  functions  of  each  liquid  are  the 
same  as  those  for  a  single  liquid  and  the  solutions  have  the  same  form. 3 

More  specifically,  the  pressure  and  velocity  perturbations  are  expressed 
as  the  sum  of  products  of  functions  of  x  and  functions  of  r.  The  functions  of 
x  are  1,  x/a,  sin  (  Xjk  x/a),  and  cos  (Ajk  x/a)  where 


Ajk  =  ( ttk/2)  [1  +  5cj]  ; 


k  =  1,3,5... (2Nk-l) 


5cj 


2) IT  (1  +  is) 


\  +  1(2-L.is)  \  Re  ■ 

(  )/3  +  is  /  1  -  i s  ;  J 


(3.2) 


(3.3) 


Re.  = 
J 


14 


(3.4) 


P,ci  =  -  [(i-s)2(x/c)(r/d)  +  :  R.  (r)  sin  (»  x/c)]  ( c/a)  K 

J  3  1  vl  K  J  * 

ujsi  *-[( i-s) (r/a)  +  (i-s)'1  r  Rjk  xk  cos  ( Ajk  x/c)]  K 
Vjsi  =  [(i-s)2  (i+s)"1(x/i)  +  Z  Ryjk  sin  (xjk  x/c)]  (c/a)  K 
wjsi  =  t"1  (i_s)2  (’+5)_1  (x/i)  +  :  «wjk  sin  (xjk  x/c)]  (c/a)  K 

V  =  Ejk  W  +  Fjk  Vr'> 

Rvjk  =  [(s-i)  a  Rjk  -  21 (a/r)  Rjk]  S_1 

R„jk  =  -  C2a  Rjk  *  Is-0  (a/r)  »Jk]  S-1 

A  A 

r  -  »jk  r/c 


S  »  s2  -  2'S  ♦  3. 


If  vj  *  u,,  the  Xjk's  for  the  same  k  are  not  the  same,  the  trigonometric 

functions  on  opposite  sides  of  the  liquid  interface  are  slightly  different, 
and  the  interface  boundary  conditions  require  a  special  least  squares  process. 
In  the  body  of  this  report,  we  will  use  a  single  ,\k  determined  by  an  average 
kinematic  viscosity,  v  . 

\k  =  (  nk/2)  [I  i  Hcl  (  ■-■a/*,1)1/'?]  (3*6) 


where 


r 


Vj  (a2  -  bx2)  +  v2  (b2  -  b2) 


a  2.2 

d  “  d  2 

In  Appendix  B,  a  much  more  accurate  treatment  of  unequal  viscosities  will  be 
described. 

IV.  THE  VISCOUS  SOLUTION 

In  Reference  3  the  viscous  parts  of  the  perturbation  variables  were  im¬ 
portant  only  in  a  small  region  near  the  cylinder  wall  and  the  endwalls,  and 
these  functions  were  computed  by  the  use  of  unsteady  boundary  layer  equations 
in  these  regions.  It  is  reasonable  to  expect  that  the  viscous  parts  can  also 
have  contributions  near  the  liquid  interface  and  that  derivatives  in  the 
radial  direction  are  much  larger  than  those  in  the  circumferential  or  axial 
directions.  Thus,  we  will  assume  that  the  viscous  functions  away  from  the 
endwalls  satisfy  the  same  equations  as  those  for  the  cylindrical  wall  boundary 
1 ayers. 3 


(s  -  t)  wjsv  *  a2  Rej-1 


2  -1  ^  u isv 

(s  - i]  uisv  - a  Rej  77- 


a  Js-V  =  ?  w 
a  9  r  1  Wjsv  • 


(4.1) 


(4.2) 


(4.3) 


The  solutions  to  Eqs.  (4.1  -  4.3)  are 


r  -  a 


-  (r  -  b, ) 


wlsv  "  w10  e 


+  W11  e 


(4.4) 


-(r  -  b,) 


w2s v  =  W21  8 


u,  =  u.n  e 
lsv  10 


+  w?2  e 


r  -  a  -(r  -  bj) 

a  6al  a  6al 

+  U11  6 


(4.5) 


(4.6) 
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nr-r.u- "■  !.■  1/ J  J 


r  -  b, 


-(r  -  b2) 


a  5 


u2sv  =  U21  e 


la2 


a  5 


Plsv  =  2  5 


'al 


+  u22  e 


r  -  a 


a2 


•(r  -  b.)  -| 


(4.7) 


a  6 


w10  e 


al 


a  6 


-  wue 


al 


(4.8) 


P2sv  3  2  6a2 


r  -  b, 


-(r  -  b,) 


a  6 


W21e 


a  2 


a  6 


where  <5 


aj 


-  w22e 


■1/2 


a2 


(4.9) 


and  the  eight  coefficients  Ujk,  Wjk  are  functions  of  x.  Four  conditions  on 

the  coefficients  come  from  the  continuity  of  tangential  velocities  and  shears 
[Eqs.  (2.17  -  2.18)]  at  the  interface. 


where 


AWsi 

+  w10  el 

+  ^  *  ^21  m  w22  e2  =  ^ 

(4.10) 

Ausi 

+  u10  el 

+  ii^2  *  ^21  "  u22  e2  =  ^ 

(4.11) 

W1Q 

el  “  W11 

=  N  (w2l  *  w22  g2^ 

(4.12) 

u10 

£1  -  U11 

=  N  (u2l  -  u22  e_2 ) 

(4.13) 

AW  . 

SI 

®-Wlsi 

-"2s,Jr  . 

r  -  b. 

AUS, 

=  Is i  ' 

'  “2S']n  -  h 

r  -  bj 

17 


*  *  .  •  J.*  '  • 

.  . 

-V 

JU 

► 

»*“. 

\ 

b1  -  a 

a  6al 

el  =  e 

_  . 

• 

b2  -  b1 

a  5a2 

e2  -  e 

R 

•/ 

„ 

N  =  P21  /  . 

1 

Four  more  conditions  come  from  the  no-slip  condition  at  the  cylindrical 
wall  [Eqs.  (2.22,  2.24)]  and  the  no-shear  condition  of  the  free  surface  [Eqs. 

(2.20  -  2.21)]  . 

1 

'« 

wlsi(a,x)  +  w1Q  +  wn  ex  =  i  (s  -  i  )K(x/a)  • 

(4.14) 

1 

* 

ulsi(a,x)  +  u10  +  un  tx  =  (s  -  i)  K- 

(4.15) 

I 

W21  e2  "  w22  =  ®  . 

(4.16) 

1 

r 

U21  e2  u22  =  0' 

(4.17) 

. 

I 

Eqs.  (4.10  -  4.17)  can  now  be  used  to  determine  the  eight  coefficient 
functions.  If  the  total  liquid  occupies  an  annular  region  that  is  thicker 
than  ten  boundary  layer  thicknesses,  the  product  e2  clu''te  small  and  can 

t 

1 

jt 

be  neglected.  The  resulting  expressions  for  the  coefficient  functions  are 
given  in  Table  2.  If  the  liquid  interface  is  more  than  ten  boundary  layer 

thicknesses  from  both  the  cylindrical  surface  and  the  free  surface,  both 
and  &2  are  ver7  smal 1  and  can  be  neglected  in  Table  2. 

'  i 

'I 

1 

FI 

Eqs.  (4.8  -  4.9)  can  be  used  in  conjunction  with  Eqs.  (4.12,4.16)  to 
yield  general  relations  for  the  viscous  perturbation  pressures  at  the  three 
boundaries. 

- 

'  1 
-■ 

« 

P2s>  (b2.*)  -  0. 

(4.18) 

1 

Plsv  (bl’x^  =  P21  p2sv  (bl**)  • 

(4.19) 

*S 

^  • 

plsv  ^a,x)  =  2  5al  (w10  "  W11  el)  • 

(4.20) 

V 
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TABLE  2.  VISCOUS  COEFFICIENT  FUNCTIONS  FOR  INTERNAL  FREE  SURFACE. 


{[1  +  N  +  (1  -  N)  e22]  ua  +  N  £l  &  u$i }  D-1 
{(1  -  N)  ejUa-  N  (1  -  e22)  A  usi}D_1 


{2  u a  +  (1  +  e^2 )  A  u$i  }  D”1 


A  u  ■  D' 


^  A  Usi 


{[1  +  N  +  (1  -  N)  e22]  va  +  N  A* }  D"1 
{-(1  -  N)  ej  va  +  N  (1  -  e22)  A*}  D'1 


{2  e, V  +  (1  +  e,2)  A*}  D_1 


A*  IT1 
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In  order  to  compute  the  viscous 

inequality  involving  the  magnitude  of 

definition  of  5  ■ 
aj 


A  P 

exp  {  -  — r  }|  =  exp 


a  6 


aj 


[-2"1/2s  r, 

1 — rrr — 


a  5 


aj 


radial  velocity,  we  must  use  an 
a  complex  exponential.  From  the 


(4.21) 


(4.22) 


where  z  =  2" I  6  .I"1  >  0 

a  1  dj  1 


Equations  (4.4  -  4.7)  can  now  be  substituted  in  the  continuity  equation. 
The  results  can  ther^  be  integrated  and  simplified  by  use  of  Inequality  (4.22) 
and  neglecting  |  j | ^  terms. 


r  'lsv  *  a6al 


r  -  a 


-(r  -  bj) 


a  6 


v10  e 


al 


a  6 


+  vlle 


al 


(4.23) 


r  v„  =  a  6  0 
2sv  a2 


r  -  b 


1 


•(r  -  b2) 


a  5a2  a  6a2 

v21e  +  v22  6 


(4.24) 


where  v1Q  =  i  «10  -  a 

V11  =  -[i  W11  -  bl 
V21  =  1  W21  *  bl  U21 

^22  ~  " ^  w22  ™  ^2  U22J 

Eqs.  (4.23  -  4.24)  can  now  be  used  in  conjunction  with  Eqs.  (4.11,  4.13  - 
4.17)  to  give  relations  for  the  radial  viscous  velocity  at  the  three  bound¬ 
aries. 


20 


V2sv  {b2*x)  =  ° 


(4.26) 


''isv  (^l*x^  =  P21  v2s v  <bl*x) 


f>21  ^  V21  +  c2  v22^  (a/bi^  5d2 


(4.26) 


vlsv  (a,x)  -  [v1Q  +  vn  £l]  5al 


(4.27) 


where  vJk  are  given  in  Table  2  (relation  (4.22)  was  used  to  simplify  entries 
in  this  table). 

V.  INV  ISC  ID  BOUNDARY  CONDITIONS 

In  the  previous  section,  we  used  eight  of  the  twelve  liquid  boundary  con¬ 
ditions  to  determine  the  viscous  perturbations  in  terms  of  the  inviscid 
perturbations.  We  will  now  completely  determine  the  inviscid  perturbation  by 
use  of  the  remaining  four  conditions — Eqs.  (2.16,  2.19,  2.23)  and  Vjs  =  v2s  at 

the  interface.  These  conditions  can  be  simply  stated  by  use  of  Eqs.  (4.18, 
4.19,  4.26,  4.26): 


at  r  -  b. 


plsi  + 


bl  vlsi 


[n  ,  bl  v2si  1 
i[p2si  a  (s  -  i )  J 


(5.1) 


vls i  =  v2si 


(5.2) 


at  r  =  bo 


at  r  -  a 


p2si  + 


b2  v2si 
a  (s  -  i 


vlsi  +  vlsv  =  f1  *  (x/a>  K 


(5.3) 


(5.4) 


eqs.  (6.1  -  6.4)  with  v^sv  and  v^sv  set  equal  to  zero  were  used  by  Scott9 

to  obtain  ei genf requencies  tor  a  completely  inviscid  liquid  payload.  As  can 
be  seen  from  Table  1,  for  each  k  four  constants  must  be  determined  to 
completely  determine  the  inviscid  perturbation  functions,  i.e.,  Ejk,  62^,  Fjk, 

•  (For  the  special  case  of  100%  filled  cylinder  (62  -  U),  is  zero  and 


> 


instead  of  4Nk  constants  only  3NK  constants  remain  to  be  determined.)  To 

obtain  these  conditions,  we  fit  (x/c)  to  a  series  in  sin  (  ^x./c)  by  least 
squares 

x/c  =  T  ak  sin  ( Akx/c)  .  (5.5) 


Next,  Eq.  (5.5)  and  the  correspondi ng  expansions  for  the  pressure  and  velocity 
functions  from  Table  1  are  substituted  in  Eqs.  (5.1-5. 4)  and  coefficients  are 
equated. 


^1  k  1  )  -  R2k  ^bl  ^  ~  ( bi / d)  t^vl k ( bl )  ”  P21  Rv2k^bl^^”^ 

(1  -  p?J  s2  (i  -  s) 

=  - — -  ( b * /a )  a 

i  +  s  ik 


(5.6) 


^vlk(bl^  "  Rv2k(bl^  “  (1  '  P21 ^ 
R2k(b2}  -  (b2/a>  <s  ‘  i)"1  Rv2k 
Rvlk(a)  •  A2  5ala  ^lk^  +  >;1 


6a2  Rk  =  ^  “  P21  ^  el  5a2  A1  ak 

<b2>  '  <Va>  *k 

N  5al  Rk*  “o'1  *  ak 


where 


R* 

< 


(d 


2a 

e9  ) 


a  -  Kz^ 


-  2  c1(a/b1)[RvU(a)  +  a  R;u(a)]}  D1 


(5.7) 


(5.8) 


(5.9) 


Rr =  2  cRvik(bi) +  bi  R;ik<bi>  ■  wbi>  •  bi  R;2k(bi^ 


r 


s  (i  -  s) 

L  1  +  s  . 


-l 


Aj  =  4  (a/bj) 


A,  -  11  *  N  +  (1  -  N )  U22  -  e^2)]  Dq_1 
D;1  =  D  -  al  [I  .  N  -  (1  -  N)(n2-  e22)] 
0  y  ~  ^  '  tl  ■*  ^  1  ( 1  t  -  ^2  )  . 


VI.  LIQUID  MOMENT 

For  coning  motion  described  by  Eq.  (2.1),  the  linear  liquid  pitch  and  yaw 
moment  is  defined  to  be:3 


"ly  +  1  MLZ  =  V2^t  (CLSM  +  iCLIM)^  e1S4> 


(6.1) 


where  m.  =  2 tt  p^a  c  . 


The  major  components  of  this  liquid  moment  are  due  to  the  pressure  on  the 
lateral  wall  and  the  endwalls  of  the  container.  Lesser  components  are  due  to 
the  viscous  wall  shear  on  the  lateral  and  endwalls.  Thus,  the  liquid  moment 
coefficient  can  be  given  as  a  sum  of  four  terms.  (For  simplicity  these  terms 
■’Ml  be  computed  for  the  center  of  the  cylinder  at  the  center  of  mass  of  the 
projectile.  More  complete  expressions  are  given  in  Reference  3.) 


c  (CLSM  +  l  CLIM)  -  +  mpe  +  myJl  +  mye  .  (6.2) 


By  use  of  Eqs.  (2.2  -  2.4)  and  Eq.  (2.10),  the  linear  periodic  part  of 
the  pressure  can  be  computed  in  cylinder-fixed  coordinates 


- —  =  R  {[p1c,  (r,x)j  es<t>  "  i0}  for  b, 

1$1  1 

=  (-'21  R  f^p2si  ^r>x^  eS'<  '  10)  for  b2 


where  pJ$1  (r,x)  =  pjsi  (r,x)  -  (rx/a2)  K. 

Eq.  (6.3)  for  the  inv’scid  pressure  and  Eq.  (4.8)  tor  the  viscous  lateral  wall 
pressure  can  be  integrated  to  yield  the  pressure  moment  coefficient  on  the 
1  a tera l  wall. 

2  3 


<  r  <  a 

(6.3) 

<  r  <-  b1 


m  =  i  (2  acK) 
P  * 


x  Lp1c,  (a,x)  +  P1cu  (a,x)J  dx 


(6.4) 


where  Plsv  (a,x)  =  2  6fll  (w1Q  -  w^). 


Since  the  viscous  pressure  on  the  endwalls  is  zero,  the  expression  for  the 
pressure  moment  coefficient  on  the  endwalls  is  slightly  simpler. 


mpe  =  "  1  (a  CK) 


Plc-i  (r»C)  r  dr 


(6.5) 


p2si  ^r,c^  r  dr 


The  viscous  moment  coefficient  on  the  lateral  wall  can  be  computed  by  use 
of  Eqs.  (4.4)  and  (4.6). 


mVJ  =  (2acK  Rej  «al)' 


Cia  (u1Q  -  ull^  +  x  ( wio  ”  wll^  dx*  (6.6) 


The  viscous  moment  coefficient  on  the  endwalls  is  a  little  more  difficult  to 
compute  since  a  change  in  kinematic  viscosity  can  occur  across  the  interface 
as  well  as  a  change  in  density.  The  relations  of  Reference  3  can,  however,  be 
used  to  obtain  the  following  result. 


v1/2 

a  K  /  1 


w^  (r)  r  dr  +  N  1  w^  (r)  r  dr 


(6.7) 


w,(r)  =  2 ( 1+i s ) ( c/a)  K  -  w.  ,(r,c)  +  i  v.  ,(r,c) 


where 


I 


>ii/;l'ssion 

:  *’  :  '  ,  t  1  ■!.  :  ,  :  4-i.  ib. 6  -  5.9)  yield  values  of  E^,  Epk, 

f  •  ,  '  v  i.  •  -  ■  :  -  '■!-  *  .  opponent  of  tne  liquid  moment.  Values 

of  s  that  mim-  tne  four  by  four  determinant  of  the  system  zero  are  called 
M  i '  ‘ 1  a  ’  ,  where  k  is  tne  axial  mode  number  and  n  is  the  radial  mode 

!i/iivr.  .re  i  riai oar  -  porn  :jr  the  s k n * s  are  frequencies,  ,  of  the 

1 1  ansi  -»■  i  r  !.•:  1  ,  ••  in  t  n,  liquid  and  the  real  parts  are  the  damping  rates  of 
that  notion.  ; for  a  single  liquid,  is  the  smallest  rkn  and  Tkn  increases 

nonotom  c al ly  with  increasing  n.)  When  the  forced  coning  motion  of  the  pro¬ 
jectile  has  a  frequency  near  an  ei genf requency ,  the  liquid  side  moment  has  a 
la.al  iiiui'i'p.  For  this  reason,  tne  eigenfrequencies  are  of  considerable 

inte''eot  to  a  project  i  le  designer.* 

The  equations  of  this  paper  have  been  coded  for  the  VAX  11/780  computer 
and  eigenfrequencies  can  be  computed  for  specified  values  of  c/a,  Re,  k,  p^., 

N,  Dj./a,  and  bp/a.  Codes  for  both  an  interior  free  surface  and  rigid  rod  are 

available.  !n  Figs.  1-3,  sample  plots  of  x^  are  given  for  a  fineness  ratio 

of  3.1  ano  Reynolds  number  of  40,000.  Each  figure  plots  x3i  versus  the  loca¬ 
tion  of  tne  liquid  interface  for  three  values  of  p^. 

In  Fiyure  1,  we  see  that  for  a  fully-filled  container,  x-^  has  the 

single-liquid  value  of  .035  at  bj/a  =  0,1  and  reaches  a  maximum  value  for  b^/a 

-  .75.  Thus,  the  maximum  change  in  ei genf requency  occurs  when  there  are 
approximately  equal  volumes  of  the  two  liquids. 

A  fill  ratio  of  98f  (bp/a  =  .141)  is  considered  for  a  free  surface  (Fig. 

2)  and  a  central  rod  (Fly.  3).  Although  there  is  very  little  change  from  the 
fully-filled  case  for  the  free  surface,  the  central  rod  has  a  strong  effect. 
For  all  density  ratios,  the  single-liquid  ei genf requency  of  .113  occurs  for 
three  values  of  bj/a  (0,  .53,  1)  and  the  maximum  effect  is  at  .8. 


•  ■  j  xv 

.  .  ,  },  r-  TO i  ' 


.  (":  -  /  ). 
i  ?:.■ !  1  .  7r,’  ' 


Scott's  experimental  data  were  for  c-a  -  3.1/.’,  Rn ,  =  ?  -  10*,  =  .82, 

and  f  =  0.9,  1.  (The  second,  third,  and  ■  >  ,i  r  *•  n  L  ,  >  1  ,t  Lett's  "able  I  are 

:ni  s- labeled  in  Reference  9.  In  his  notation,  "n>\,  s'n/.i  he  identified  as 

-2  -2 

d/c,  bdc  ,  and  adc  .  )  His  data  are  compare-1  w'"-  t"e  tneory  of  this  paper 
in  Figure  A,  and  we  see  that  the  maximum  error  i  *  1 .  I*'  Reference  3  it  was 
found  that  much  better  fits  could  be  obtained  by  a  slightly  different 

value  of  the  fineness  ratio.  Indeed,  effective  ;alit-s  of  fineness  ratio  which 
were  0.5%  greater  than  the  measured  value  would  give  excellent  agreement.  In 
Figure  4  it  is  shown  that  a  0.7%  greater  fineness  ratio  (c/a  =  3.150)  gives 
excellent  agreement  for  six  different  experiments . 


An  important  characteristic  of  an  ei genf reguency  is  the  occurrence  of  a 
maximum  side  moment  coefficient  for  a  coning  motion  with  constant  damping  and 
a  frequency  near  the  ei  genf  requency .  The  theory  of  this  report  was  used  to 
compute  the  maximum  side  moment  coefficients  for  the  conditions  of  Figure  1 
and  constant  amplitude  coning  motion  with  frequency  near  t^.  The  resulting 

curves  are  given  in  Figure  5  and  have  a  number  of  interesting  properties. 


As  b^/a  varies  from  0  to  1,  the  heavier  liquid  is  replaced  by  the  lighter 
liquid.  Since  the  side  moment  is  defined  in  terms  of  the  density  of  the 
heavier  liquid,  a  requirement  for  our  C^sm  (bi/a)  i s  that 

cl$mU)  =  o21  CLSM(0)  •  (7*D 


This  requirement  was  satisfied  by  our  calculations.  It  is  important  to  note 
that  a  large  change  in  side  moment  occurs  when  the  liquid  interface  is  near 
the  cylindrical  wall  and  the  interface  boundary  layer  overlaps  the  wall 
boundary  layer.  Therefore,  small  amounts  of  heavier  liquid  can  have  a  large 
effect  on  the  side  moment  exerted  by  a  spinning  liquid  payload. 


An  interesting  feature  of  the  =  .4  side  moment  coefficient  curve  is 

the  presence  of  two  small  peaks.  These  are  caused  by  the  equality  of  and 

t52  f°r  these  l°cations  of  the  interface.  Figure  C  compares  the  complete 

dependence  of  x^  on  interface  location  with  that  of  x^  for  p^^  =  .4.  For 

n  =  1  the  ei  genf  requency  has  one  maximum,  while  for  n  =  2  the  corresponding 
ei genf requency  has  two  maxima, and  for  small  enough  density  ratio  those  eigen- 
frequency  curves  can  intersect. 


The  coalescence  of  ei genf requencies  for  special  interface  locations  also 
occurs  for  the  fineness  ratio  of  4.29,  which  has  been  extensively  studied  by 
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S' 


rc 


Me  i  .  .-rc  ,  >r  i  ji,.  eie  in  ,ping  ''it.'  prolned  by  a  mixture  of  oil  and 
nit"  ;  uie  Mr.vr  .  1  •  Tor  a  mixture  of  159  silicon 

■ :  :  e  •  ■  >•  •  •  •  p’ott-.-.i  versus  in  Figure  111.  This 

••i  *•  ,  •  ;  -  ,i  ,  i  ■  ir  i  ■:  c.,-:t.ii  ier  mounted  in  the  gyroscope.  The 

.t '  -.,i  u-  pi  i.  Hi  .id!  v  e  .;-.e  ;  to  predict  the  damping  rate,  and  this 

,.i r  1 1  ti  /  is  sho.vn  ai  t. he  j a sned  curve.  The  shape  of  the  curve  is  good  but 

t  u.’i-o  i  •,  a  frequency,  Mas  or  .9025.  This  Mas  corresponds  to  changing  the 
fir..-', ess  mi. ii>  0.3-  from  the  actual  v  1 1  ue  of  3.126  to  an  effective  value  of 


Till.  SUMMARY 

•  '■  ’  r-  i  Hy  layer  t  neor ,  uas  ieci  developed  for  a  two-liquid  payload 

i:  .  it  i  •,  .  ..;i(  •  j  an  i  spinning. 

i 

Tr.e  pi.",  r.e  :  e  l  u:o  f  reque  ;.i  '■  es  and  side  moments  agree  well  with  available 
t.-r’ me  a..-,  e-'iept  fur  an  unexplained  frequency  bias. 

i.ar jt;  .  nar-ies  in  si  le  mo-ients  can  he  caused  by  the  presence  of  a  small 
dunum  ,r  oaa  ■  ;er  1  igu  i  i. 

» 

:•)<  .  a  oral  i  ■  >  n  c.  of  tne  liquid  interface,  pairs  of  eigenfrequencies 

a1  ■  ;  'a  1  •.  ••  . 
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and  Various  Density  Ratios. 
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Figure  7.  t7?  versus  b./a  for  Re,  =  10°,  c/a  =  4.29,  f  =  1  and  Various  Density  Ratios 


Figure  8.  t  and  xn  _  versus  bj/a  for  Rej  =  106,  c/a  -  4.29, 
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1 1- h LC  r  hi-  CENTRAL  ROD 


t.rfn.r  of  cfntkal  rod 


N  1  l  . 


The  effect  ..r  .1  <  en •  r  j i  red  op  .1  single  liquid  payload  has  been  discussed 


1  n  del  ji  1  by  !-f  >/:  er .  •' 
conditions  si-m  !.ji  1  0 


.  n  1  i ) on  ,  , 2.1)  -  2.21)  at  r  =  b2  are  replaced  by 
-  •' .  24  '■ . 


! s  -  i )  (Op/a)  K 


v 2 s ^ b 2 ’ x )  "  ~l-s  '  (x/d)  K  • 


w./  (b.,,x)  -  i(s  -  i)  (x/a)  K. 


Equation  (A2),  uf  course,  replaces  tq.  (5.3)  as  the  boundary  condition  at  r  = 
Do.  Eqs.  (4.1b  -  4.1/)  are  then  replaced  by 

Wpsi(b2,x)  +■  w21e2  +  w22  =  i  (s-i)  K(x/a)  (A4) 

U2si(b2’x)  +  U21£2  +  U22  =  *(b2/a)  ‘  (A5) 

A  revised  version  of  fable  2  can  be  computed  for  the  viscous  coefficient 
functions  (see  Table  Al).  For  these  viscous  coefficient  functions  the  only 
radial  viscous  velocity  to  be  affected  is  that  at  r  =  bp  [Eq.  (4.25)]. 


2sv'  2 


(b2,x)  -  (a/bp)  ( vpp  +  c2  Vp,)  &  2  . 


The  inviscid  boundary  conditions  now  differ  from  Eqs.  (5.1 
condition  at  r  =  b?.  They  are 


-  5.4)  by  the 


at  r  -  bi 


,  bl  v lsi  f  bI  v2si  1 

1 1  s  i  a" 7 s  '  -  T)  ■‘■'21  ^2s i  a  (s  -  i )  J 


■  r  <■  :!  v.  * .  Scott,  "Py-iamias  of  a  Liyuii-Filled  Shell: 
\  i  r  •/  Rod,"  Ballistic  Research  Laboratory, 

.r->n>:  1,  "aryl  i»;  /,  RRL  Memo  ran  iwn  Report  3RL-VR- LtRl , 


;  ■.!  }  ti  i-Fi!  led  Shell:  Viscous  Effects  t»;  u 
Ro.i hill  isti-c  Research  Laboratory, 
f .  in  :  ,  BRL  '-U’M orandun  Report  BRL-PR-ldSO, 
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TABLE  Al.  VISCOUS  COEFFICIENT  FUNCTIONS  FOR  CENTRAL  ROD 


w1()  =  {[1+N  -  (1-N)  ^  wd  +  H  A  wsi}D”1 

=  {(1-N)  -  N(1+e22)a  w$.+  2  N  £2wb}U 

W21  =  ^  elwa+  U  +  e]^)A  wSi"  ^(^“N)  wb  ^  ^ 
w22  =  ^-e2  A  wsi  +  +  vl'N)  Ej2]  w^}  D  ^ 

u1Q  =  {[1+N  -  (1-N)  c22]  ua  +  N  £l  a  u$i }  D'1 
u 1 1  =  {(1-N)  e^u a~  N  (l  +  c22)  A  us-+2  N  e2  ub^  * 
U21  =  ^  elua+  ^  +  ei^  A  usi"  e2  ub)D  1 

u22  =  ^“e2  A  usi+  ^-1+N  +  ^1-N)  el?^  ub^_1 

v1Q  =  {[1+N  -  (1-N)  e22]  va  +  N  EjA*}!)"1 

vn  =  (-(1-N)  C;lvd+  N(l+e22)  a*  -  2  N  e2vb}D-1 

V21  =  ^2  elVa  +  ^  +  Gl^  A*  "  e2^~N)  vb^”^ 
v22  =  (e2  A*  -  [1+N  +  (1-N)  e^]  vb ID”1 


where 

D  =  1  +  N  +  (1  -  N)  (ej2  -  e22) 
wb  =  i  (s  -  i)  (x/a)  K  -  w2$i(b2,x) 
ub  =  (s  -  i)  (b2/a)  K  -  u2si(b2/x) 
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1  tit?  difference  between  tqs.  (All)  and  (A13)  is  only  important  when  b£  is  very 

bit, ail  while  the  difference  between  Eqs.  (A12)  and  (A14)  is  unimportant  for 
Reynolds  numbers  appropriate  to  boundary  layer  theory. 

Arter  the  j -i  jt=  1  substitutions  of  the  Fourier  series  expansions,  new  forms 
of  lqs.  (o.o  -  b.9)  tor  the  tjk's  and  Fjk's  can  be  obtained. 
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^ Ik ( b  1 )  *  P21  R2k(bl5  '  (bl/a'  (s  '  LHvU(b1)  -  p21  Rv2k(bi)] 


(1  -  p21)  (i  -  s) 

m 


( b 1 /a )  ak 


Rvlk(bl^  ‘  Rv2k(bl)  ‘  ^  “  P21  ^  6a2  Rk  =  ^  "  p21^el  ”  e2^  6a2  Alak 


Rv2k^b2^  +  a  °a2  D  ^2  Rv2k^b2^  +  (a/b2^  e2  6a2  ^2  1  Rk* 


-  -2$-U  -s)a, 
i  +  s  k 


R,lk<a>  -  a  sal  A2  "ilk  (■»  *  5al  N  “o'1  Rk’  '  --TXk"?)  ak  • 


rfhere 


“f  ■  <a  <>  *  ci2  *  e22>  [R;ikibi>  -  R;2k<bi>j 


-  2  (a/bj)  ex  [Rylk(a)  +  a  R^lk(a)] 


+  2  (a/bj)  e2  CRv2k(b2)  +  b2  R;2k(b2)]}  D'¬ 


'S  (i  -  s ) i  n  -1 


\  • 4  <a/V  [Hir1r1]  V 


A2  -  [1  +  N  -  (1  -  N)(£l2  +  e22)]  D0_1 


Dn  =  D  -  6,,  [1+N  -  (1-N)(e  2+  e,2)] 


Dj  =  D  -  5a2  (l-p21)(l+ei2+  e22)(a/b1) 


n2  =  D  f  dd2  D  (a/b2) . 

The  linear  liquid  moment  coefficient  for  a  cylinder  with  a  central  rod 
can  be  expressed  as  the  sum  of  six  terms — the  four  terms  of  Eq.  (6.2)  plus  the 
pressure  and  viscous  contributions  from  the  central  rod. 


i  ( C.  Clul  +  i  C.  m)  =  m  +  m  „  +  m  +  m  +  in  +  m 
L3M  LIM'  p  £  pe  pr  v£  ve  vr 


( A19 ) 


The  two  additional  terms  are  quite  similar  to  mpJi  and  myJi,  respectively. 


"'pr  =  -i  b2  p21  (2a2cK) 


c 

"L 


x  ^p2si^b2,x)  +  p2sv^b2*x^  dx  , 


( A20) 


where  p2$v  -  -2  <5a2  (w22  -  w21  e2\ 


(2  a  cK 


Ci  b2  ^22*  e2u 2 1  ^  x  (w22"e2  ^21^  dx  ^ ^ ^ ) 


In  Figure  A1  the  maximum  side  moment  coefficient  associated  with  is 

plotted  versus  b^/a  for  the  central  rod  of  Figure  3.  Note  the  sharp  changes 

in  side  moment  for  the  interface  near  the  central  rod  and  near  the  outer 
cylindrical  wall.  The  two  bumps  in  the  p2^  =  .4  curve  are  caused  by  the 

coalescence  of  the  and  the  t^2  eigenfrequencies. 
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APPENDIX  B.  EFFECT  OF  DIFFERENT  KINEMATIC  VISCOSITIES 


As  was  stated  in  the  body  of  this  report,  the  approximation  give"  in  Eg. 
(3.5)  can  be  avoided  by  a  much  better  least  squares  approximation.  Exact 
expressions  for  PjS-j.  vjsi>  dnd  vjsv  can  he  drived  from  the  series  of  Table  1 


and  the  other  relations  given  in  the  main  part  of  this  report,  and  these  are 
tabulated  in  Table  Bl.  Next,  these  expressions  for  the  pressure  and  radial 
velocity  perturbation  functions  are  substituted  in  Eqs.  (5. 1-5.4)  to  obtain 
four  equations  for  the  E k ^ ‘ s  and  Fk|'s: 


kj 


l  [Elk  sik  *  Flk  “ft3  sin  Ulk  x/c)  *  E  [E2k  c;k  *  F2k  11  iJ  5,n  ('2k  </c> 


'  d;  */c 


£  =  1,2, 3, 4 


(Bl) 


k  =  1 ,3,5,7 , . . ,2Nk  -  1 


where  the  coefficients  A^,  8^,  C^,  0^,  d^  are  given  in  Table  B2. 


For  Nk  sets  of  the  parameters  (E^,  E2k ,  Fn<»  ) ,  these  four  equations 

can  not  be  satisfied  exactly.  We  can,  however,  consider  the  squared  sum  of 
the  residuals  for  each  equation  produced  by  a  particular  selection: 


V  j  _c  I  E  ^Elk  Akk  +  Flk  B£k^  sin  ^ikx/c^ 


+  ^  ^2k  ^ £k  *  ^2k  ^ik^  x/c)  -  d0  x/c  dx 


2k 


(B2) 


We  seek  those  values  of  the  parameters  Egk,  E2k,  Flk,  F2k  that  minimize 

(in  a  least  squares  sense)  the  R^s.  For  a  single  liquid  and  100%  fill  (b2  = 

0),  the  Fk '  s  satisfy  the  condition  at  r  =  b2  while  the  Ek's  satisfy  the 

condition  at  r  =  a.  For  two  liquids,  the  conditions  at  b2  and  a  are 

represented  in  Eqs.  (81-B2)  by  t  =  3  and  4,  respectively.  Hence,  in  analogy 
with  the  simpler  case,  we  select  the  F^'s  to  minimize  Rg  and  the  Ejk's  to 

minimize  .  To  satisfy  the  interface  conditions  (i  -  1,2),  we  select  the 

Fjk's  to  minimize  R^  and  the  E2k's  to  minimize  R?.  The  resulting  4NK 

equations  are  given  in  Table  B3. 


Tables  Bl  and  B2  have  to  be  modified  for  a  fully-filled  cylinder  with  a 
central  rod  of  radius  b2.  The  pressure  and  inviscid  radial  velocity  functions 
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!■>  7 3 tile  til  are  unchanged,  but  three  viscous  radial  velocity  functions  have  to 
'■v  derived  from  Table  A1  and  are  given  in  Table  B4. 


A,  was 

indicated  in 

Appendix  A,  the  boundary  conditions  at  r  =  for 

the 

surf  ace 
ent 1 r  1  e 

(Eq.  (5.3'' 

1  by  -  3. 

is  replaced  by  Eq.  (AD.  The  corresponding  Eq. 
Tne  new  coefficients  are: 

(Bl) 

A3, 

'a  3 

h? 

Ik 

’’  k 

'  'd/b2'  V 

h2k 

'  Ik 

=  f?K  (V  f 

( a/h? )  sa2  h3k 

(B4 ) 

i)3k 

=  y2k  (h2)  + 

(a/b2)  5a2  h*k 

l13 

--  2  s ( i  -  s) 

(i  +  s)'1  -  (a/b?)  3a2  h*  . 

lainimj  coefficients  for  it,  =  1,2,4  are  given  in  Table  B3  when  the  hp. 


nn< 1 s  if  Table  B4  are  used. 


As  a  comparison  of  approximation  (Eq.  3.5)  with  the  better  approximation 
of  tois  appendix,  the  maximum  side  moment  coefficient  for  T33,  Re  =  40,000, 

c/a  =  3.1,  by  =  0,  ,)  =  0.6  is  plotted  versus  bj/a  in  Figure  Bl.  Comparison  of 

tne  two  approximations  are  given  for  kinematic  viscosity  ratios  (v^  = 

of  0.1,  1,  and  10.  We  see  that  approximation  (Eq.  3.5)  underestimates  the 
effect  of  unequal  kinematic  viscosities. 


TABLE  Bl.  PRESSURE  AND  RADIAL  VELOCITY  FUNCTIONS. 
(i-s)2(x/c)(r/a)  +  Z  [EjkJlk(f)  +  F.J^tr)]  sin  (\jkx/c)}  (c/a)  K 
-s)2  ( i +s) _1 ( x/c)  +  Z  [Ejkfjk(r)  +  Fjkgjk(r)]  sin  (xjkx/c)}  (c/a) 

=  {h ^ x/c  +  L  [E^k  b^k  +  E^k  ^2^ ]  sin  (x^kx/c) 

'  [E2k  h3k  +  F2k  h4k3  Sin  (X2kx/c)}  (c/a)  6al  K 
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J  dr 
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J  dr 
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TABLE  B4 .  VISCOUS  RADIAL  VELOCITY  FUNCTIONS  FOR  CENTRAL  ROD. 


,(a,x)  =  Lhj  x/c  +  Z  [Elk  hlk  +  Flk  h2k]  sin  (Alkx/c) 
+  '  ^E2k  h3k  +  F2k  h4k^  sin  (A2kx/c'}  6al  * 


2sv(bl’x)  =  {h2x/c  +  E  [Elk  hbk  +  Flk  h6k]  Sin  (Alkx/c) 


+  ^E2k  h/k  +  F2k  h8k-*  sin  ^x2kx/c^}  (c/bi^  6a2  K 


2svv  2 


(b?,x)  =  {hf  x/c  +  Z  [Elk  h*k  +  Flk  h*k]  sin  (A^x/c) 


+  ^E2k  h3k  +  F2k  h4k-^  s1n  (A2kx/c^}  (c/b2^  6a2  K 
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•2  t,  l)"1  (a)  +  (I  +  e,  +  k  ‘)  D"1  fb, 


lsu)ijb  RADIAL  VELOCirv  FUNCTIONS  FOk  CENTRAL  ROD  (Continued) 


9 1 i:  ( <J )  f  !  ’  •  f  '  i?c)  o"1  y lk ( ^  1 ) 


A  n-l 


*■  r2  1  u  f'z<{l'l)  +  1  *1  11  1  f~2i^nz: 


■l  *  •/)  *  2  *2  I'*1  92k(»2> 


ir1  f1K(t>l) 


,:'1  9lK(nl> 


2  4  (1  -  “>  ‘  I  -  2  ‘2  'WV 


2  ,/  11  -  N)  U*1  J  ij  (b  )  -  2  ,  l)'1  .],  (»,) 


I 

■-j 


APPENDIX  C 


AX  ISYMMt  TR  IC  EIGENVALUES 


PACE 


APPENDIX  C.  AX  SYMMETRIC  EIGENVALUES 
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made 
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Re  t e  i 
The  n 


wh  ef  i 

the 
i  n  1 1  >  i 


1  h i  e i ijeriv.) I aes  or  tnis  report  were  determined  by  linear  perturbation 
.P-,  whose  azimuthal  variation  was  exp  (-io).  A  more  general  expression 
.  i-i  -i  j  since  the  boundary  conditions  have  an  azimuthal  variation 
s f  r :  ;  t. .  i  -i  -  1,  the  steady-state  response  required  a  consideration  of 
-.a  I  ue  of  .;i  wnly.  Experimental  measurements  of  eiyenfrequencies  have  been 

f.  =  r  j.h  v.un-t  ric  waves  (m  -  D)  in  a  single  liquid.1-1  In  this  appendix  we 
o,o  ,  ohi ■  <;,ju at  ions  for  axi symmetric  waves  in  two  liquids. 

[he  mviscid  perturbation  functions  for  any  m  have  been  derived  in 
>rii.e  t:y.  In  Table  Cl  the  appropriate  expressions  are  given  for  m  *  1. 
lore  general  definitions  for  -S  and  6  .,also  given  in  Reference  C2,  are 

J  Cl 


w  1  f  J - Hi:  *l/--  (Cl) 

1  1,1  /2(m+is)  J 


=  -L E? ( ii»+  is)]  1  [  ( zf -in- is)  oj|n  *  -  (2+m+is)  3jm  ^  1  (C2) 

-  (c/a)  [s  -  (2+m)  i]1/2  ke.1/2  (C3) 

-  (c/a)  [s  +  ( 2-in)  iji/2  Re.1/2,  (C4) 


;  the  ctmiplex  roots  are  selected  to  have  positive  real  parts. 

Gince  velocities  at  rigid  boundaries  are  zero  for  in  t  1  perturbations, 
corresponding  viscous  coefficient  functions  can  be  computed  for  an 
dial  free  surface  from  Table  2  when  wa ,  u^,  and  va  are  defined  to  be 


-  j  s  ,  ( a  ,  x ) . 


(C5) 


■i‘i float  i.n:  of  the  Inertia  Osoillatior.s  of  a 
HRL  Contract  Report  2?St  September 
lo-'ph’jr,.  Aiit.rophja.  Fluid  Pynami oe ,  Vol. 
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TABLE  Cl.  INV ISC  ID  PERTURBATION  FUNCTIONS  FOR  m  *  1. 
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Table  A1  can  be  used  to  compute  viscous  coefficient  functions  for  a  central 
rod  and  m  t  1  if,  in  addition  to  Eqs.  (C5-C7),  similar  definitions  are 
assigned  to  w^,  and  v^: 


wb  =  -  w2si  (b2,x) 
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r  a  b2  • 


ueneia!  versions  at  Kqs.  (b.b-b.y)  for  in  v  1  can  now  he  derived 
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f  i  1  i  vUnr  1'  v2Knr  1'  1  v2kinv  1 ' J  ‘ 

t  p. .  ‘  li-  i  4  ■  inive  been  coded  for  m  =  0  and  eigenvalues  can  be  computed. 
4l  ! 1  id  '  •  t  .  i.  i  ->d  rn  ei  jenf  requency  measurement  for  Re  =  43,000 

an  t  i  •»  ,  |n  M  |, I,-  '  !  ,,  ver-.us  bj/d  is  plotted  for  ^  -  .4,  .6, 

.  .  -  i  •  - . I  r  il  l.  f  I  single  liquid  IS  plotted  for  bj/d  -  0,  1. 


Re  =  43000,  c/a  =  0.995, 


Figure  Cl.  "^lO  versus  for  Re ^  -  4.3  x  10^,  c/a  =  0.995,  f  1,  m  0  and 

Various  Density  Ratios. 
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radius  of  the  cylindrical  cavity  containing  the  liquids 

coefficients  in  a  least  squares  fit  of  x/c  to  a  series  in 
sin  (  \  x/c)  ,  Eq.  (5.5) 

radius  of  the  interface  of  the  two  liquids  for  K}  =  0 

radius  of  an  air  core  (or  central  rod)  for  =  0 
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pressure  moment  coefficient  on  the  end  walls,  Eq.  (6.5) 

pressure  moment  coefficient  on  the  lateral  wall,  Eq.  (6.4) 
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eigenvalues  of  s:  the  values  of  s  that  make  the 
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time 


coefficients  in  the  expansion  of  Ujsv,  Eqs.  (4. 6-4. 7) 
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axial  coordinate  in  the  inertial  system 

axial  coordinate  in  the  aerobal 1 i Stic  system 
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Bessel  function  of  the  first  kind  of  order  1 
i Jt(0) 
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amplitude  or  the  coning  motion 
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1 

'2 

''J 


V 
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( K i /K i )  ,  non-dimensional  i zed  damping 
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exp  [(b2-b1)/a  <5a2] 
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